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Abstract 

Non-minimally coupled scalar field models of dark energy are equiv- 
alent to an interacting quintessence in the Einstein's frame. Consider- 
ing two special important choices of the potential of the scalar field, i.e. 
nearly flat and thawing potential, one has an analytical expression for 
the equation of state parameter as a function of the density parameter 
of the scalar field for any choice. Here we investigate the non-minimal 
quintessence model by applying the method of statefinder diagnosis 
to it and plotting the evolutionary trajectories of the statefinder pa- 
rameters. 



1 Introduction 

According to the supernova observations [TJ, our universe has an accelerating 
expansion at present epoch [2j [3]. Accepting general relativity theory as 
the underlying theory for describing the dynamics of the universe, there is a 
missing energy-momentum component of the universe, dubbed "dark energy" 
[I], providing the present acceleration. There are several candidate for such 
a dark energy, cosmological constant, chaplyian gas, quintessence and so on 
(See Ref jl] and references therein). 

Some cosmological models involve an interaction between dark matter and 
dark energy [5] . This removes the need for a fine-tuned cosmological constant 
to get the ratio of dark matter to dark energy of order unity nowadays. This 
is widely known as coincidence problem [6]. 

A particular class of these models is Brans-Dicke (BD) type scalar field 
model. One can deal with BD scalar field model in two frames," Jordan 
frame" and "Einstein frame". In Jordan frame, the the scalar field does not 
appear in the action of matter fields but couples non-minimally to gravity. It 
is possible to make a particular conformal transformation, converting the the- 
ory to a scalar-tensor model in which the scalar field couples conformally to 
matter and minimally to gravity. Describing the scalar field as quintessence, 
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hereafter we shall refer to this model as conformally coupled quintessence 
model (CCQ). There has been a lot of discussion about the choice of one of 
these frames as the physical frame [TJ [8] . Here we pay our attention to the 
Einstein frame. This frame is a suitable framework to make the model con- 
sistent with the solar system constraint and further observational constraints 
from Big-Bang nucleosynthesis and inflation [H] . 

Following Caldwell and Linder [10] , the scalar field models of quintessence 
can be divided into two categories, called "freezing" and "thawing" models. 
In the former, the equation of state parameter has a decreasing behavior but 
in the latter, it has a value near —1 initially, and then increases with time to 
less negative values. 

In [TTJ, it is shown that non-interacting thawing quintessence model 
with nearly flat potentials (the potentials satisfying the slow-roll conditions: 



(f c^J ^ 1' VlyF ^ 1) P rov ides a natural way to produce a value of cu near 
— 1 today. Generalizing this to CCQ model with nearly flat potentials [T2] , 
shows that there exists a universal behavior for uo which is different from 
thawing behavior initially. In a recent paper [13] . some conditions on the 
potential of the scalar field are derived which is different from the slow-roll 
conditions initially and lead to the thawing behavior for all times. 
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On the other hand as demonstrated in Refs. [2] , it is possible to discrim- 
inate different models of dark energy from each other by some parameters, 
called the statefinder parameters (r, s) , firstly proposed in [15] and [16] . The 
first of these, r, is the jerk parameter and the other is a function of jerk and 
the decelerating parameters. The statefinder pair depends on the metric of 
space-time and is constructed using the second and third derivatives of the 
scale factor as: 

T — 1 

S= 3(g-l/2) (2) 
where q = — a /a 2 is the deceleration parameter. 

In this paper we will discuss further the nearly flat and thawing CCQ 
models of dark energy by statefinder diagnostic. In the next section, we 
will start with BD interacting dark energy model. According to [12] and 
|13j . application of some conditions on the potential, nearly flat and thawing 
conditions, leads to an analytical expression for u in any case. Then, we 
apply the statefinder diagnostic to CCQ model in section 3. In the last 
section we will give some conclusions. Throughout this work we have chosen 
the units 8ttG — c — 1. 
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2 The model 

The general action of BD scalar tensor theory in the Jordan frame is: 

Sj = J d 4 x^g~ - - 2C/($) + . (3) 

where R is the Ricci scalar of the metric g^, U{$) is the potential of the 
scalar field, uj is the BD coupling constant and C m is the matter Lagrangian. 
Under the conformal transformation g^ v = e^g^ , in which In $ = and 
( = 3^;, one arrives at the action of BD theory in the Einstein frame: 

S E = J <?xy/=g[R - l -(V V f - V{<p) + C m (e-^ 9flu )}. (4) 

where V(<p) = e~ 2 ^U{^)). 

Now, consider a spatially flat FRW universe occupied by the pressureless 
matter. The cosmological equations of motion are: 

H 2 = -(p^ + p m ) (5) 

H = -^\Pm + P v + P v ] (6) 

and the scalar field evolution is governed by the following equation of motion: 

(p + 3H (p + V v = s^pp m (7) 

where p m , p^ = \(p 2 + V{ip) and p v = |<^ 2 — V((p) denote the energy density 
of cosmic fluid, the energy density and pressure density of the scalar field in 



the Einstein frame respectively and (3 = y §C From equation (j7j), one can 
easily see that the energy density of the scalar field satisfies the following 
conservation law: 

p\ + 3#(1 + u)p v = \J^pipp m (8) 

in which u = p v / p v is the equation of state parameter for the scalar field. 
Combining the above equations, the continuity equation for the cosmic fluid 
can be derived as: 

p m + 3Hp m = -^pLpp m (9) 
By integrating equation (EJ), one obtains: 

/ \ / a \ 3 (v-yo) , . 

Pm(t) = POm — e (10) 

Veto/ 

in which p 0m , a and (p are the current values of the matter density, the scale 
factor and the scalar field respectively. Therefore the usual dependence of the 
non-relativistic matter upon the scale factor is modified by an exponential 
factor due to the interaction with the scalar field. 

Taking the time derivative of V((p) = ^-^- P? anc ^ us i n g the continuity 
equation (jSJ), one obtains [TT] : 



J = -3(1 -uj 2 ) 



_v^ + i i-% 



3(l+a;) V V VQ ^ 



where Vt m and Q v are the density parameters of matter and the scalar field 
respectively (Qi = pi/3H 2 ) and prime denotes derivative with respect to 
In a.. Ignoring the interaction of the scalar field and the cosmic fluid, this 
equation reduces to the corresponding equation of Steinhardt et al. [TS] for 
a non interacting scalar field. 

Following the same method as refs [HI [T2], using the relation <fi 2 = (1 + 
o>)/ty, the continuity equation (jSJ) reads as [13]: 



which is a useful relation to derive the redshift dependence of the cosmological 
quantities in the next section. 

Moreover by using equation (flQj) . one can express equations (JSJ) and (J6]) 
in the following form, convenient for the construction method proposed by 



where x = 1 + z in which z = a /a — 1 is the redshift parameter and a is the 
present scale factor of the universe. These are two coupled equations allowing 
one to reconstruct the potential V(z) and the scalar field (p(z) knowing fi 0m , 
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(12) 



[19j 




(13) 




(14) 



ifo and also using the observed H(z) from the luminosity distance. 

Since the left hand side of equation (j!4p is non-negative, the dynamical 
expansion of the universe is restricted by the following inequality: 

^ > 3H 2 n 0m (l + z)*e-K*-*> (15) 

This is the weak energy condition for BD cosmology in the Einstein's frame. 
We see that although it is not possible to express the Hubble parameter as 
an explicit function of the cosmological redshift in this model, there is a 
necessary condition on the H(z) in this case. Ignoring the interaction and 
thus the exponential factor, inequality f|T5|) reduces to what is derived by 
Sahni and Starobinsky in j2] for non interacting case. 



3 Statefinder diagnostic of CCQ model of dark 
energy 

According to the previous section, using BD theory of gravity in the Ein- 
stein's frame leads to a scalar-tensor theory in which there is a particular 
interaction between the non-relativistic matter and the scalar field. It is a 
well known fact that the coupled scalar field model with — 1 < uj < — 1/3 
(called usually interacting quintessence) has the capability of explaining the 
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current cosmic acceleration [20] • However to make these models acceptable, 
some limitations on the form of V has to be set, depending on the form of 
the interaction term [T8| |2T| [22]. 

In JT3] , the necessary conditions for the existence of thawing behavior for 
CCQ model is found. Assuming thawing behavior, the following relation, 
called thawing condition, should be satisfied 



which shows that it is necessary that A increases with time when (p and 
Q v are increasing functions. Dividing equation ffTTT) by equation ffT2l . one 
arrives at a differential equation for u as a function of [13]. Replacing 
A with expression (fT6|) and retaining terms to the lowest order in 1 + w (w 
is near —1), the differential equation of uo is exactly solvable. The resulting 
analytical expression expression for the state parameter is as follows |13j : 



2A 

I I — w . \ 

1 +UJ 



1 _ O \ " n 9)i 1 / 2 + j4 i o 



(17) 

where A is a positive constant, xo is an integration constant, A = 1+/3^J^\ 
and 2 -Pi is the Gauss Hypergeometric function. 

In [13] it is shown that in the case of CCQ model, A neither is a small 



value nor is a constant. Therefore the nearly flat potentials do not lead to the 
thawing behavior. However recently, study of the behavior of the equation 
of state parameter for nearly flat potentials has attracted a lot of attention 
[TTltrS]. The authors of [12] have shown that the equation of state parameter 
firstly increases with time and then approaches asymptotically to a value 
near to —1. As mentioned in [12], assuming the slow-roll conditions for the 
potential, one can show that y < 1 which ensures that A is approximately 



constant up to now, i.e. 

A^A = - — 

where Aq is a small constant evaluated at the initial value of <pq. Combining 



(18) 



equations ( fill) and ( Tl2|) and making two assumptions, the first one is that uo 
is near —1 and the second is that the condition ( {TBI is satisfied, these yield 
again to a differential equation for w as a function of Q v [12] which gives the 
following analytical expression for the equation of state parameter: 



1 +u = [_£_ - ( -I - i) H*> - ^ In l±Xi | - a]f (19) 



"If 

where a = — ^ 2 ^^ n ' — 2y/ ^^ Qi in which f2j is some small initial value of fl v 
such that oji = —1. 

From equations §5§ and (jSJ), it is straightforward to show that the decel- 
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eration parameter takes the following form as: 



q=-(l + 3uQ v ) (20) 



After differentiating equation (jSJ), using relations (jSJ) and (J9]), one finds: 

9 3 3 / — 

r = 1 + -uj(1 + cj)^ - -u/^ - —=ujyjl + ujj%{\ - Sl v ) (21) 
Z Z ZyZ 

Inserting ([2D]) and fl2D in (fJD gives: 

s = l + w-- ^_ y vTT^ 22 

3^ 3^ ^ 

which explicitly depends on in contrast to the non-interacting quintessence. 

Let us now have a detailed look at how the statefinder parameters behave 
for thawing and nearly flat CCQ model of dark energy. As we have seen 
before, in these cases, it is possible to derive an analytical expression for the 
statefinder pairs as a function of the density parameter of the scalar field, Q^, 
without considering a special form for the potential. To do this, one can use 
the relation (fTTI) . the thawing condition (fISl) and the analytical expression of 
u>, equation (IT7|) . for the case of thawing CCQ model and the corresponding 
relations (TTT|) . (fTSl) and (IT91 in the case of nearly flat CCQ model. 

The time evolution of the statefinder pairs (r, s) for thawing CCQ model 
has been shown in Figure [T] with Ao = 0.9. This value of Ao is chosen such 
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that u has a value near —1 today. Also the constant Xo is chosen such that 
the initial condition uj = — 1 holds for = 0.001 and /3 has been set equal 
to 0.5. KCDM (A-cold dark matter) universe corresponds to the fixed point 
(1,0). We see that the evolution of trajectories of statefinders pairs pass 
from the point (r, s) ~ (1, —0.07) in the past when z ~ 4.62 and 0,$ ~ 0.01. 
And then after passing A CDM fixed point, r decreases whereas s increases 
to the point (r, s) ~ (0.04, 0.27) at z ~ -0.74 (f2 ~ 0.99) in the future. The 
location of today's point is (r,s) ~ (0.57,0.15) when ~ 0.7. This shows 
the present 'distance' of thawing CCQ model of dark energy from KCDM 
model. The time evolution of the pairs (r, q) is indicated by Figure [2j The 
dashed, thick and thin curves have A = 0.8, 0.9, 1 respectively. We see that 
both A CDM and thawing CCQ models start evolving from the same point, 
(r, q) = (1, 0.5) which corresponds to SCDM (standard cold dark matter) 
universe. In A CDM scenario, the evolution is along a horizontal line ends 
at SS (steady-state) fixed point, (r, q) = (1,-1), corresponding to de Sitter 
expansion. However in our model, q has a decreasing behavior whereas the 
value of r first increases and then monotonically decreases. The trajectory 
goes to (r, q) ~ (0.04, —0.68) at z ~ —0.7 in the future. 

In Figures [3] and HI we have shown the deceleration parameter and the 
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equation of state parameter as functions of redshift parameter. These figures 
have been plotted numerically using the relations ffl2l . (|T7|) and (1201) . 

In Figure [5] we have shown time evolution of the statefinder pairs for 
nearly flat CCQ model of dark energy with Ao = 0.4. This small value of 
Ao ensures that the variation of the potential during the evolution of the 
universe is very small. Assuming that for an initial value of f^, say fij, 
Qi = 0.001, we have to = —1. This determines the constant a in the relation 
f fT9|) setting (3 = 0.5 here. The (r, s) trajectory has two branches. We find 
that the first branch comes asymptotically from r ~ 1, s — > — oo goes to 
r ~ 1.33, s — > +oo for a change of in the interval [0.001,0.17]. Another 
branch comes along r ~ 1.33 asymptote, passes from A CDM point and then 
goes to (r,s) ~ (0.78,0.05) at z ~ -0.7 ~ 0.99) in the future. The 
divergent behavior and discontinuity of s occurs at ~ 0.17 and is due to 
vanishing u; at this point. Along the first branch, both r and s are increasing, 
however for the other branch r monotonically decreases whereas s increases. 
The location of today's point is (r, s) ~ (0.94, 0.02) when ~ 0.7. This 
shows that this model has less distance from A CDM model in comparison 
to thawing model. This result is satisfied in both models independent of the 
allowed chosen value of A . 



13 



The evolutionary track in (r, q) plane is shown in Figure O The tick and 
thin curves are corresponding to Ao = 0.4 and 0.1 respectively. It is started 
from SCDM point as the same as the A CDM evolutionary path, but first 
r decreases and q increases slightly, then after passing a period in which r 
and q have increasing behavior, those decreases monotonically. This period 
corresponds to very high redshift. The trajectory goes to (r, q) ~ (0.78, —0.9) 
at z ~ —0.71 (Qfj, ~ 0.99) in the future. At the end, we present the plot of q 
and ui with respect z. 

4 conclusion 

In this paper, we have applied the statefinder diagnostic to the CCQ model 
of dark energy. And we have plotted the trajectories in the (r, s), (r,q), 
(q, z) and (a;, z) for nearly flat and thawing potentials. As it is apparent 
from Figures [3] and the deceleration parameter decreases monotonically 
with redshift for both the nearly flat and thawing potentials, however it goes 
to more negative values at future in the case of nearly flat model. Moreover 
nearly flat potentials force the equation of state parameter changes in a wide 
range with redshift whereas it remains near —1 for thawing model, as one 
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Figure 1: Plot of r versus s for thawing model. 

expected. This is exactly what one finds from Figures H] and El Also the 
value of r decreases with redshift steadily for both kinds of potentials as well 
as non- interacting quintessence model [15]. 
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Figure 2: Plot of r versus q for thawing model. 
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Figure 3: Plot of q versus redshift for thawing model. 
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Figure 4: Plot of u versus redshift for thawing model. 




Figure 5: Plot of r versus s for nearly-flat model. 
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Figure 6: Plot of r versus q for nearly-flat model. 
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Figure 7: Plot of q versus redshift for nearly-flat model. 




Figure 8: Plot of oo versus redshift for nearly- flat model. 



18 



References 

[1] A. G. Riess, et al., Astron. J. 116, 1009, (1998) S. Perlmutter, et al, 
Astron. J. 517, 565, (1999) 

[2] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9, 373, (2000); 
V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 15, 2105, (2006) 

[3] T. Padmanabhan, Phys. Rept. 380, 235, (2003) 

[4] E. J. Copeland, M. Sami, S. Tsujikawa, Int. J. Mod. Phys. D, 15, 1753, 
(2006); P. J. E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559, (2003) 

[5] L. Amendola, Phys. Rev. D 62, 043511 (2000);W. Zimdahl, D. Pavon 
and L. P. Chimento, Phys. Lett. B 521, 133,(2001); W. Zimdahl, Int. J. 
Mod. Phys. D 14, 2319 (2005); D. Pavon and W. Zimdahl, Phys. Lett. 
B 628, 206 (2005); L. P. Chimento, A. S. Jakubi, D. Pavon and W. 
Zimdahl, Phys. Rev. D 67, 083513 (2003) 

[6] P. J. Steinhardt, in Critical Problems in Physics, edited by V. L. Fitch 
and D. R. Marlow, Prinston Univ. Press, Prinston, NJ, (1997) 

[7] S. Capozziello, S. Nojiri, S. D. Odintsov, A. Troisi, Phys. Lett. B 639, 
135143 (2006) 

19 



[8] V. Gorini, A. Kamenshchik, U. Moschella, V. Pasquier, A. Starobinsky, 
Phys. Rev. D 72, 103518 (2005) 

[9] J. A. R. Cembranos, Phys. Rev. D 73, 064029 (2006); A. A. starobinsky, 
JETP Lett. 86, 157 (2007); T. Fulkner, M. Tegmark, E. F. Bunn and 
M. Mao, Phys. Rev. D76, 063505 (2007) 

[10] R. R. Caldwell, E. V. Linder, Phys. Rev. Lett. 95, 141301, (2005) 

[11] R. J. Scherrer and A.A. Sen, Phy. Rev. D, 78,083515, (2008) 

[12] A. A. Sen, G. Gupta and S. Das, arXive:0901.0173 (2009) 

[13] M. Roshan and F. Shojai, Phys. Rev. D 80, 043508 (2009) 

[14] V. Sahni, A. Shafieloo and A. A. Starobinsky, Phys. Rev. D 78, 103502 
(2008); C. Zunckel and C. Clarkson, Phys. Rev. Lett. 101, 181301 
(2008); C. Feng, Phys. Lett. B 670, 231 (2008); D. J. Liu and W. Z. 
Liu, Phys. Rev. D 77, 027301 (2008); H. Wei and R. G. Cai, Phys. Lett. 
B 655, 1 (2007); Z. L. Yi and T. J. Zhang, Phys. Rev. D 75, 083515 
(2007); M. R. Setare, J. Zhang and X. Zhang, JCAP 0703, 007 (2007); 
M. G. Hu and X. H. Meng, Phys. Lett. B 635, 186 (2006) 



20 



[15] V. Sahni, T.D. Saini, A. A. Starobinsky, U. Alam, JETP Lett., 77, 201 
(2003) 

[16] U. Alam, V. Sahni, T. D. Saini, A. A. Starobinsky, MNRAS, 344, 1057 
(2003) 

[17] M. Roshan and F. Shojai, Phys. Rev. D 79, 103510 (2009) 

[18] P. J. Steinhardt, L. Wang, and I. Zlatev, Phys.Rev.Lett. 82, 896 (1999); 
P. J. Steinhardt, L. Wang, I. Zlatev, Phys.Rev. D 59, 123504 (1999) 

[19] A.A. Starobinsky, JETP Lett. 68,757 (1998), T. D. Saini, S. Raychaud- 
hury, V. Sahni, A. A. Starobinsky, Phys.Rev.Lett. 85, 1162 (2000) 

[20] D. Wands, E. S. Copeland and A. Liddle, Ann. N. Y. Acad. Sci. 688, 647 
(1993); L. Amendola, Phys. Rev. D62 043511 (2000); J. M. F. Maia and 
J. A. S. Lima, Phys. Rev. D 65, 083513 (2002); W. Zimdahl, D. Pavon 
and L. P. Chimento, Phys. Lett. B 521, 133 (2001); L. P. Chimento, A. 
S. Jakubi, D. Pavon and W. Zimdahl, Phys. Rev. D 67, 083513 (2003); 

[21] T. Chiba, Phys. Rev. D 66, 063514 (2002) 

[22] R. Das, T. W. Kephart, R. J. Scherrer , Phys. Rev. D 74, 103515 (2006) 



21 



